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An outstanding problem in the field of high-transition-temperature (high Tc) supercon¬ 
ductivity is the identification of the normal state out of which superconductivity emerges 
in the mysterious underdoped regime.^ The normal state uncomplicated by thermal fluc¬ 
tuations is effectively accessed by the use of applied magnetic fields sufficiently strong to 
suppress long-range superconductivity at low temperatures.^’^ Proposals in which the nor¬ 
mal ground state is characterised by small Fermi surface pockets that exist in the absence 
of symmetry breaking^’have been superseded by models based on the existence of a su¬ 
perlattice that breaks the translational symmetry of the underlying lattice.’"^^ Recently, a 
charge superlattice model that positions a small electron-like Fermi pocket in the vicinity 
of the nodes (where the superconducting gap is minimum)*’^’has been proposed a re¬ 
placement for the prevalent superlattice models^®"^'* that position the Fermi pocket in the 
vicinity of the pseudogap at the antinodes (where the superconducting gap is maximum).^* 
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Although some ingredients of symmetry breaking have been recently revealed by crys¬ 
tallographic studies, their relevance to the electronic structure remains unresolved. 
Here we report angle-resolved quantum oscillation measurements in the underdoped cop¬ 
per oxide YBa2Cu306+x- These measurements reveal a normal ground state comprising 
electron-like Fermi surface pockets located in the vicinity of the superconducting gap min¬ 
ima (or nodes), and further point to an underlying superlattice structure of low frequency 
and long wavelength with features in common with the charge order identified recently by 
complementary spectroscopic techniques.^'*’ 

The normal ground state eleetronie strueture revealed by our measurements is summarised in 
Fig. 1. Quantum oseillations measured in the eontaetless eleetrieal resistivity of the underdoped 
euprate YBa2Cu306.56 are shown in Fig. 2a. The prominent oseillatory beat strueture reveals a 
sizeable frequeney spread.^’ Measurements are made as a funetion of magnetie field up to 85 
T over a wide range of field orientations as defined by the polar angle 6 and the azimuthal angle 
(j). This data may be used to identify eharaeteristies of the momentum-spaee eleetronie strueture 
in the way that, for example, X-ray data may be used to identify the real-spaee lattiee strueture 
of a erystal. In partieular, it ean be employed to identify the eorreet Fermi surfaee model. 
In the ease of underdoped YBa2Cu306+x, whieh has a primitive orthorhombie lattiee. Fig. 2b 
illustrates the quasi-two dimensional eleetronie strueture eharaeterised by a eylindrieal Fermi 
surfaee with fundamental neek and belly warping that would be expeeted in the absenee of a 
superlattiee (Fig. 3a-d). A distinguishing eharaeteristie of this neek and belly Fermi surfaee is 
that it gives rise to a resonanee in the amplitude of quantum oseillations at a polar angle known 
as the Yamaji angle.^^ The resonanee in amplitude is expeeted to grow with the size of the 
frequency spread, and its location is expected at approximately is expected at approximately 
60° in YBa 2 Cu 306 +x given the measured mean diameter of the cylinder corresponding to the 
dominant quantum oscillation frequency, and the c-axis lattice constant (Fig. 2b,c). Strikingly, 


2 



however, despite the sizeable frequeney spread of the quantum oseillation speetrum (Fig. 2a), 
higher angle data in Fig. 2a reveal the Yamaji resonanee antieipated for a fundamental neek and 
belly Fermi surfaee geometry to be absent, pointing to an alternative Fermi surfaee geometry. 
Supporting analyses are shown in Methods and Extended Data Figs. 1,2,3,4. 

To arrive at the eorreet Fermi surfaee model we are guided by the following key experi¬ 
mental findings: (1) as shown in Fig. 4a, the replaeement of the expeeted Yamaji resonanee 
indeed by an anti-resonanee; (2) as shown in Fig. 4b, the fourfold anisotropy in 0 dependenee 
of the quantum oseillation amplitude for different values of Q, with a maximisation of quantum 
oseillation amplitude along the a and b erystallographie direetions and minimisation along the 
diagonal direetions, for example, 0 = 45°; and (3) as shown in Fig. 4b; the enhaneement in 
fourfold anisotropy as a funetion of azimuthal angle with inereasing values of polar angle 9. 

While the observed experimental features are in sharp varianee with those expeeted for a 
dominant fundamental neek and belly Fermi surfaee geometry, we find that a dominant stag¬ 
gered twofold Fermi surfaee with warping maxima along the diagonal direetions in momentum 
spaee (Fig. la-c. Fig. 3f-h) prediets the experimental features that we have observed. The sim¬ 
ulation of the quantum oseillation waveform (Fig. 3h, and Fig. 4, and Methods) shows that sueh 
a staggered twofold Fermi surfaee geometry yields first an anti-resonance in the quantum os¬ 
cillation amplitude at a special value of 9 along the diagonal directions in 0, second, a fourfold 
anisotropy in the quantum oscillation amplitude as a function of <p that is minimised along the 
diagonal directions, and finally an enhancement in the fourfold anisotropy as a function of 0 
with increased values of 9. Layered materials with examples of both fundamental neck and 
belly and staggered twofold Fermi surface geometries are discussed in Methods. 

Simulations of a staggered twofold Fermi surface geometry with modulation amplitude 
AFtwofoid ~ 15 T, agree well with the measured quantum oscillation data as a function of 
all three parameters B, 9, and 0 (Fig. 4, Extended Data Table 2). In addition to the ampli- 


3 



tude damping factor for a staggered twofold geometry eonventional damping faetors 

are ineluded that are used to deseribe layered materials, whieh arise from thermal smearing, 
impurity seattering or quenehed inhomogeneities, Zeeman splitting and magnetie breakdown. 
The quantum oseillatory frequeney spread is eaptured by Fermi surfaee splitting from a finite 
bilayer or spin orbit eoupling. More details of the simulation are given in Methods, Extended 
Data Figs. 2c,4,5,6,7, and Extended Data Tables 1,2. 

The staggered twofold Fermi surfaee geometry we observe would not be expeeted to dom¬ 
inate within the primitive orthorhombie unit eell of underdoped YBa2Cu306+x- Yet just sueh a 
geometry would arise from the unique symmetry of the eomer T point in the Brillouin zone of a 
body-eentred orthorhombie unit eell defined in Figs. 1 and 3e-h (see Methods). At this speeial 
comer T point, the twofold in-plane rotational symmetry alternates by 90° between adjaeent 
symmetry planes. The diagonal orientation of the maximal warping direetions of the Fermi sur¬ 
faee (shown in Figs, la-c and identified from the azimuthal anisotropy in Fig. 4b), reveals the 
new body-eentred orthorhombie Brillouin zone to be oriented eoneentrieally with the original 
primitive orthorhombie Brillouin zone, and to be defined by orthogonal ordering waveveetors 

and Q 2 . The staggered twofold Fermi surfaee poekets at the eorner T point are therefore lo- 
eated in the nodal region of the original Brillouin zone. A nodal Fermi surfaee poeket eontained 
within the reeonstrueted Brillouin zone is not only eonsistent with the present measurements, 
but also with the low observed value of linear speeifie heat eoeffieient in high magnetie fields^ 
and strong ehemieal potential oseillations inferred from previous quantum oseillation studies. 

We now eonsider the possible origin of the superlattiee responsible for the above nodal 
staggered twofold Fermi surfaee poekets, whieh emerge from a reeonstruetion of a large Fermi 
surfaee determined from band strueture, eharaeteristie of the original primitive orthorhombie 
lattiee in the normal state.^ Numerous proposals for superlattiees in the eopper oxides have 
been put forward, for example in refs.^’’“^^’^"^^^^ Of partieular relevanee to our observations is 
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the charge order recently detected in YBa2Cu306+3; by techniques such as X-ray diffraction/^’^*’ 
ultrasound/’ nuclear magnetic resonance/"* and optical reflectometry/^ These observations 
point to a superlattice characterised by the same ordering wave vectors Qi and Q 2 identified by 
our measurements and defined in Figs. 1, with the superlattice scaling parameters (5i ^ ^2 ~ 0.3, 
and anisotropic amplitudes in some cases. We note that for Fermi surface reconstruction, the 
superlattice need not be strictly long-range or static, but it must not be fluctuating over a range 
much smaller than the cyclotron radius, nor with a frequency much larger than the cyclotron 
frequency. 

Numerical calculations indeed show that this type of charge order can give rise to nodal 
Fermi surface pockets similar to those observed in our measurements, potentially accompanied 
by antinodal gaps in the electronic excitation spectrum^’(schematic shown in Extended Data 
Figs. 6, 7). Furthermore, the resulting nodal electron-like Fermi surface pockets yield the ob¬ 
served negative sign of Hall coefficient in the accessed high magnetic field limit (see Methods, 
refs.^’**) in contrast to the positive value expected for the large unreconstructed Fermi surface 
in the absence of the superlattice. Our measurements do not distinguish between charge order 
with s-wave symmetry and unconventional symmetry such as d-wave symmetry (i.e., bond or¬ 
der). Further, they are compatible with charge modulation components (of wavevectors Qi and 
Q 2 ) not only of similar amplitudes, but also of significantly different amplitudes, as in a nearly 
uniaxial structure.^’ 

The strength of magnetic fields in which the present measurements were carried out (up to 
85 T in Fig. 2a and up to 100 T in Extended Data Eig. 3b) is adequate to suppress long-range su¬ 
perconducting order, and thus reveal the normal ground state underpinning superconductivity in 
underdoped YBa2Cu306+a;. The resulting normal ground state is a Eermi liquid characterised by 
staggered twofold Eermi pockets that emerge naturally from a long-range static or short-range 
slowly fluctuating superlattice defined by ordering wavevectors Qi and Q 2 (Eig. 1), support for 
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which is found from recent complementary experiments.^"^’Superconductivity may there¬ 
fore be viewed as emerging from the pairing of quasiparticles on Fermi pockets we locate in 
the nodal region of momentum space. Our finding clarifies observations from complementary 
experiments, for example Raman spectroscopy,^*’ which show that in the underdoped regime, 
Bogoliubov quasiparticles are confined to momentum space islands around the nodal regions of 
the Brillouin zone. This is in sharp contrast with the emergence of Bogoliubov quasiparticles 
from both nodal and antinodal regions of the starting large unreconstructed Fermi surface in the 
overdoped regime. 

Method Summary 

Quantum oscillations in the electrical resistivity were measured using a contactless technique 
on a high quality de-twinned single crystal of YBa 2 Cu 306.56 over a wide range of polar and 
azimuthal angles, Q and 0, in magnetic fields up to ~ 85 T and at 1.5 K (Fig. 22 and Extended 
Data Fig. 1). Measurements up to 100 T at 6^ = 0 are presented in Extended Data Eig. 3. 

The model used for the quantum oscillation simulations is discussed in the Methods. Of 
particular importance here is the Eermi surface geometrical damping or warping factor, taken 
to be of the form = Jq sin 2(j)J2{kF{c'/ 2 ) tan 6^)], which is expected for 

a staggered twofold Eermi surface geometry in a body-centred weakly-orthorhombic unit cell 
(Eig. la,b). Here kp ^ 0.13 A~^ (corresponding to the dominant frequency Fq ^ 530 T) is 
the average radius of the Eermi surface pockets in the basal plane, AFtwofow is a measure of 
the magnitude of the staggered twofold warping and Jq and J 2 are the zeroth and second-order 
Bessel functions, respectively. The results of the simulation are shown and explained in Eig. 4 
and in Extended Data Table 1. 

Eor a Eermi surface with fundamental neck and belly warping geometry in a primitive 
weakly-orthorhombic unit cell (see Eig. 2b), a different warping damping factor /?°eck-beiiy _ 
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Jo(/i:Fctan 6*) is required. The results of a simulation using with 


model parameters listed in Extended Data Table 3 are shown in Fig. 2c and Extended Data 
Fig. 2b and c. The striking difference in the angular-dependence of the quantum oscillation 
waveform corresponding to /?°eck-beiiy (illustrated in Fig. 3d and h respectively) is 

shown in Extended Data Fig. 2. 
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Figure 1: Fermi surface of YBa2Cu306 56 inferred from quantum oscillation measure¬ 
ments. a, Fermi surfaee poekets in the nodal planes within the reeonstrueted body-eentred 
orthorhombie Brillouin zone (in red outlines) and the original primitive orthorhombie Brillouin 
zone (in blue outlines), b, Staeked view of reeonstrueted Brillouin zones defined by ordering 
waveveetors Qi = 27r(±^,0, ±^) and Q 2 = 27r(0, ±^, ±^). c, Cut of b through = 0; 
the dotted line represents the Fermi surfaee eross-seetion at a eut through = tt/c. d, Blue 
surfaees represent the inferred sinusoidal superlattiee that is staggered between bilayers: here 
the unprimed and primed lattiee parameters relate to the real-spaee unit eell and the reeon¬ 
strueted superlattiee unit eell respeetively. (5i and 82 are given in the text, and the hole doping 
for YBa2Cu306.56 is approximately 10% holes per Cu atom. 
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Figure 2: Quantum oscillations defining the Fermi surface of YBa2Cu306.56* a, Oscillations 
of the scaled contactless electrical resistivity as a function of for different magnetic field 
(B) orientations defined by 0 and 6 as indicated. The upper inset defines the angles 9 and 0 
relative to the crystallographic axes, b, Schematic showing the degeneracy in the cyclotron 
orbit cross-sectional area that yields an amplitude enhancement at the Yamaji angle (green) for 
a neck and belly Fermi surface geometry (corrugation is accentuated throughout for clarity), c, 
Expected Yamaji resonance in the vicinity of ~ 60° were YBa2Cu306.56 to be described by a 
neck and belly Fermi surface geometry (see Methods, Extended Data Eig. 2b, c). Data in a and 
simulations in c have been scaled by exp for visual clarity. 
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Figure 3: Fermi surface and geometry-dependent quantum oscillation damping for differ¬ 
ent crystal structures, a-d, Primitive orthorhombie strueture. e-h, Body-eentered orthorhom- 
bie strueture. Shown are the real spaee unit eells (a, e), euts through the Brillouin zone showing 
loeal symmetries about the vertieal lines (b, f), loeation of the warped Fermi surfaee within the 
Brillouin zone in whieh symmetry points F, Z, M, and T are indieated (c, g), and the angular 
dependenees of the eorresponding quantum oseillation amplitude damping faetor R„ (d, h). 
The amplitude for the fundamental warping (d) exhibits a maximum versus 6 — known as the 
Yamaji resonanee — in eontrast to an anti-resonance for the staggered twofold warping (h). 
has been simulated using parameters in Extended Data Table 1 and a representative value of 
i? = 62 T. For the definition of and further analyses including the spin factor see the 
Methods. 
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Figure 4: Quantum oscillation data compared with a staggered twofold Fermi surface 

model, a, Real eomponent of the eross-eorrelation between the quantum oseillation data 
over a fixed range of B cos 6 for a range of measured 6 angles with a simple sinusoid 
cos{2'kF/B cosO 4- 0). F and 0 are matehed to the periodieity and phase of the oseillations 
at 6* = 38°, where a single frequeney dominates the measured quantum oseillations. The dotted 
lines show the expeeted eonstant magnitude of the maximum and minimum amplitude for an 
ideal two-dimensional Fermi surfaee. A suppression in amplitude (anti-resonanee) instead of 
a Yamaji resonanee is observed as a function of 6 (see Methods, Extended Data Fig. 2). b. 
Fourfold anisotropy in the 0-dependent amplitude of quantum oscillations for 6 = 45° and 60°. 
c. Measured quantum oscillations shown in Fig. 2 (black lines). The magenta lines in a and 
b show a simulation of and the coloured lines in c show a simulation of 4/twofold 

respectively for a staggered twofold Fermi surface model (for a value of AFtwofoid = 15 T). For 
definitions of the staggered twofold damping factor the spin damping factor R^, the 

quantum oscillation function 4/twofold, the modulation amplitude of staggered twofold geometry 
AFtwofoid, and parameter values used for the simulation, see the Methods. 
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Methods 


Experimental Details. 

High quality detwinned single crystals of YBa2Cu306.56 of dimensions 0.5 x 0.8 x 0.1 mm^ 
used in this study were grown and prepared at the University of British Columbia.Quantum 
oscillations were measured at the National High Magnetic Field Laboratory (Los Alamos and 
Tallahassee) using the contactless resistivity technique (described elsewhere^^’^^). The sample 
and the proximity detector oscillator coil to which it is coupled were rotated in situ through 
different polar {6) and azimuthal (0) angles before each magnetic field pulse extending to 65 
or 85 T. Quantum oscillation measurements for 6^ = 0 were extended to 100 T, in the regime 
where long range superconductivity is destroyed at low temperatures.^’Additional az¬ 
imuthal angle dependences were measured using a two-axis goniometer in the 45 T DC hybrid 
magnet. In all cases the sample temperature was maintained close to T ^ 1.5 K throughout the 
experiment by direct immersion in superfluid "^He. A worm drive rotator powered by a stepper 
motor at room temperatures was used for sample rotation, with a secondary angular calibration 
provided by a pancake projection coil wound in the plane of the sample. For the pulsed-field 
data, the measured component of the applied magnetic induction B cos 6 (where B ^ h-qH) 
projected along the crystalline c-axis shown in Extended Data Fig. 1 yields an uncertainty of 
approximately 0.2° or less in the sample orientation. 

The polar angles accessed in the rotation study shown in Figs. 2c and 4a and Extended Data 
Eigs. 1 and 2 are 0 = 0°, 1.3°, 11.3°, 12°, 16.3°, 18°, 21.3°, 26.3°, 31.3°, 36.3°, 38°, 41.3°, 
45.2°, 46.3°, 48°, 49°, 49.4°, 50.1°, 50.6°, 51.4°, 51.5°, 52°, 52.3°, 52.5°, 52.9°, 53.1°, 54.4°, 
54.9°, 55.5°, 56°, 56.2°, -56.95°, 57.2°, -57.4°, -58.15°, 58.2°, -59.4°, 59.6°, 60.6°, 61.2°, - 
61.4°, 61.7°, 62.5°, 62.6°, -62.7°, -63.2°, 63.4°, 63.7°, -64.1°, 64.5°, 65.5°, 66°, 66.3°, 68.1°, 
69.4° and 70.6°. Negative 9 angles refer to a measured equivalent (180-|6'|) angle, as shown in 
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Extended Data Fig. 1. 


Absence of experimental Yamaji effect in underdoped YBa2Cu306+x- 


Extended Data Fig. 2a shows a compilation of quantum oscillations measured for YBa2Cu306,56 
(also shown in Figs. 2a and 4c). In previous experiments over a limited range of magnetic field 
and orientation (indicated by thick grey dashed line in Extended Data Fig. 2b), experimental 
data was compared with a simulated quantum oscillation waveform from two cylinders with 
fundamental neck and belly geometry: 


'hneck—belly — E dj [RsRbRtR w] j 


Sin 


5 cos 6* 




( 1 ) 


where j = 1, 2 for two cylinders, and the relative amplitudes aj, phases 7 j, quantum oscilla¬ 
tion frequencies Fj, and fundamental neck and belly warpings AFneck-beiiyj were considered 
as independent parameters.The spin damping factor Rg, Dingle damping factor Rd, thermal 
damping factor Rt, and the geometrical damping, or warping factor are defined in the fol¬ 
lowing sections. In this case, the quantum oscillatory frequency spread is modelled by neck and 
belly warping, which predicts a dramatic enhancement in quantum oscillation amplitude at a 
special Yamaji polar angle of approximately 60 o. This Yamaji resonance is a consequence of 
all the orbits becoming degenerate at this special angle.^^’^^ Simulations made using parameters 
from ref.^^ (shown in Extended Data Table 3) are shown in Fig. 2b and c and Extended Data 
Fig. 2b and c. In the present experiment in underdoped YBa2Cu306+x, where an extended mag¬ 
netic field and polar and azimuthal angular range are accessed, despite the sizeable frequency 
spread, we unexpectedly find the absence of any quantum oscillation amplitude enhancement 
in the vicinity of a polar Yamaji angle of approximately 60°, with instead a suppression of the 
quantum oscillation amplitude at a special polar angle in the vicinity of ~ 60°. 
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Simulations using staggered twofold Fermi surface geometry. 

Fig. 4 shows a simulation of a Fermi surface with staggered twofold geometry that agrees with 
the experimental quantum oscillation waveform and amplitude remarkably well over the en¬ 
tire experimental range in B, 6 and 0. The simulation uses a quantum oscillation function of 
the form expected for a bilayer-split nodal Fermi surface with staggered twofold warping (see 
Fig. la-c) 

6 

twofold ^ ^ [-^w-^MB-^s-^D-^t] j COS 

J=1 

Here, oq is the amplitude prefactor (which is taken to be the same for all orbits), i?MB is the 
magnetic breakdown amplitude reduction factor to be defined in the following sections, and 
Nj counts the number of instances the same orbit is repeated within the magnetic breakdown 
network. In our model, Rd, and Rt are taken as the same for all orbits. 

Conventional damping parameters included in simulations. 

Quantum oscillation simulations include conventional thermal, Dingle, and spin damping fac¬ 
tors of the same form used for previous comparisons with quantum oscillations measured in 
the underdoped copper oxides and other layered families of niaterials.^’^’^°’^®’^’’^^’^^’^^^^^ The 
thermal damping factor is given by 


Bcos9 


— n 


( 2 ) 


Rt = ^ 

sinh Xj 

(where Xj = 2n‘^k-Bm*gjT/heB), and the Dingle damping factor is given by 

Rb = exp r "), 

\ B cos 9 ) 

(where Aj is a damping factor).The spin damping factor is given by 
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where cos 9 for a given orbit ‘j’ is determined by the projection B cos 6 oi B per¬ 

pendicular to the planes (that is, the projection parallel to the c-axis in YBa2Cu306+a;). The 
anisotropic effective g-factor has the form g^- = which can be renor¬ 

malised with respect to the electron g-factor by spin orbit coupling,many body effects^^ and 
the presence of small band gaps.^^’^^ Here, m*^ and g^ refer to the respective values of m*gj and 
gl, when B is parallel to the crystalline c-axis, while is the anisotropy in the spin 

susceptibility. Because of the multiple frequencies in the model, it is not possible to identify the 
g-factors uniquely. The values of g* and here represent parameters used for the simulation. 

The quantum oscillatory frequency spread is modelled by a splitting of the Fermi surface, 
which can arise not only by tunnelling between bilayers,but also, for example, via spin- 
orbit effects under certain conditions.This leads to two starting frequencies that are denoted 
for reasons that will become clear below, as Fq — 2Fgpiit and Fq -f 2Fspiit.^*’ Magnetic breakdown 
tunneling (in the nodal region where the splitting is smallest) can then give rise to a series of 
combination frequencies. Here we consider cyclotron orbits corresponding to a nodal bilayer- 
split Fermi surface from charge ordering shown in Extended Data Fig. 6 and 7 (modelled in 
refs.,^’^’^^’^°’^^ ). The combination frequencies from the Fermi surface orbits depicted in Ex¬ 
tended Data Eig. 7 are listed in Extended Data Table 2. Eor a small magnetic breakdown gap, a 
high tunnelling probability P causes the orbit of frequency Fq to dominate the quantum oscilla¬ 
tion amplitude (as seen in Extended Data Eig. 3a). The beat structure in the quantum oscillation 
simulation is caused primarily by the superposition of the dominant Fq oscillations and weaker 
amplitude oscillations of frequency Fq — AFgput and Fq -f AFgput. Oscillations of frequency 
Fo — 2 AFgpiit and Fq -f 2 AFgpnt are expected to be even more strongly attenuated. The magnetic 
breakdown amplitude reduction factor has the form 

fmb = {tVpy‘'{VT^y\ 
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in which and count the number of magnetie breakdown tunnelling and Bragg refleetion 
events en route around the orbit, having transmitted amplitudes iy/P and yjl — P respeetively. 
The magnetie breakdown probability is given by P = exp(—Po/-Bcos6'), where Pq is the 
eharacteristie magnetic breakdown field. 

Angle-dependent damping term from Fermi surface warping geometry. 

The Fermi surface warping geometry in layered materials leads to a quantum oseillation damp¬ 
ing faetor with an azimuthal (0) and polar (6*) angular variation that is very sensitive to the 
energy band dispersion perpendieular to the planes (known as the ‘warping’ of the Fermi sur- 
faee). The amplitude damping faetor (P„) whieh arises as a eonsequenee of weak Fermi surfaee 
warping is obtained by (1) an expansion of the Fermi surfaee waveveetors in terms of cylindrical 
harmonics, and retaining only the lowest order harmonic in kz, (2) an evaluation of the Fermi 
surfaee cross-sectional area A{kz,9,(l)) in a plane normal to the field orientation defined by 
{9, 0) and erossing kz, and (3) via a kz integration of exp{ihA{kz, 9, 0)/(eP)). 

For the case of a primitive weakly-orthorhombic unit cell^^ (with lattice dimension c, see 
Fig. 3a) in the limit of weak warping, we thus arrive at the quantum oseillation damping faetor 
for leading order neek and belly geometry, in the simplified form 

^neck-belly ^ J tan 0)1 (3) 

P COS0 

In eontrast, for the ease of a body-eentred weakly-orthorhombic unit cell (with lattiee dimension 
c' = 2c, see Fig. Id and 3e) in the limit of weak warping, the quantum oseillation damping faetor 
for leading order staggered twofold geometry has the simplified form 

P“ = Jo [^^^sin 20 J 2 (A:i.(c 72 )tan 0 )l (4) 

P cos o 

The former has well known examples oeeurring in materials families ineluding organie con- 
ductors^^^^ while the latter oeeurs at the Brillouin zone eomer in materials families ineluding 
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body-centered tetragonal pnictides,^^ ruthenates^^ and Tl-based (overdoped) high T^. supercon- 
ductors."^^’^^ 

Here the effective Fermi wavevector k-p (the average radius of the Fermi pocket in the basal 
plane) and warping parameters AFn are related to the cylindrical harmonic expansion parameter 
k^i, defined in ref.^^ by kp = /cqo and = /i/i;oo^ 2 n,i/e, as is appropriate to the leading order, 
z/ = 1, cylindrical harmonic in kz. In this work, we refer to the amplitude of fundamental neck 
and belly geometry AFq as AFneck-beiiy, and the amplitude of staggered twofold geometry 
AFi as AFtwofoid- The relevant Fermi surface geometry depends on the local symmetry in 
the Brillouin zone. Examples of an isotropic azimuthal anisotropy for a fundamental neck and 
belly Fermi surface geometry and a fourfold azimuthal anisotropy for a staggered twofold Fermi 
surface geometry are shown in Fig. 3. 

While a fundamental neck and belly geometry would be expected to dominate at all locations 
in a primitive orthorhombic Brillouin zone, a staggered twofold Fermi surface geometry would 
arise from the unique symmetry of the special corner T point in the Brillouin zone of a body — 
centred orthorhombic unit cell (shown in Figs. 3f,g, and la-c). At this comer T point, the 
twofold in-plane rotational symmetry is seen to be rotated by 90° between adjacent symmetry 
planes. The special corner point is characterised by a leading order expansion of the Fermi 
wavenumber 

kp{kz, 0) = kp{(t)) + A/ctwofoid cos{{c!/2)kz) sin(20), (5) 

in cylindrical harmonics, where kz is the wavenumber along the c-axis, c' = 2c (Fig. Id and 3e), 
and the parameter A/ctwofoid is the amplitude of modulation along the c-axis, with the following 
symmetries: (1) invariance under the joint transformation kz ^ kz + 27r/c' and 0 —)■ 0 -f 90°, 
(2) mirror symmetry about fcz = 0 and (3) mirror symmetry about 0 = 45° (that is, the nodal 
planes in YBa2Cu306+x shown in Fig. la). While this special symmetry of the corner point 
does not support a neck and belly geometry, it instead supports a Fermi surface with a stag- 
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gered twofold geometry (Figs, la to c and 3g, and examples from other layered unconventional 
superconductors^^’ ^ i • 52^ 

We see from the simulated geometrical quantum oscillation amplitude damping factor as 
a function of polar and azimuthal angle that a fundamental neck and belly geometry and a 
staggered twofold Fermi surface geometry would exhibit strikingly different angular depen¬ 
dences (shown in Figs. 3a-d and 3e-h respectively, simulated using equations 3 and 4). A stag¬ 
gered twofold Fermi surface geometry would yield an amplitude suppression (anti-resonance) 
in quantum oscillation amplitude at a special polar angle, in contrast to the enhancement in 
quantum oscillation amplitude expected at a special Yamaji polar angle^^’^^ for a fundamen¬ 
tal neck and belly Fermi surface geometry. Furthermore, while a fundamental neck and belly 
Fermi surface would yield an isotropic azimuthal dependence of quantum oscillation amplitude 
to leading order, a staggered twofold Fermi surface geometry would yield a fourfold azimuthal 
dependence of quantum oscillation amplitude to leading order. 

Staggered twofold Fermi surface geometry simulation. 

Extended Data Table 1 shows the value of parameters used for the quantum oscillation wave¬ 
form simulated for a staggered twofold Fermi surface using Eqns. 2 and 4, which is compared 
with the experimental data in Eig. 4 and Extended Data Eigs. 2c and 4. Eor simplicity, the 
value of magnetic breakdown field Bq, warping AFt^ofoid, damping A and effective mass is 
taken to be the same for all orbits (enabling us to drop the subscript j). Eurthermore, only two 
sets of anisotropic g-factors are considered: orbits Fi, F 2 , F 4 , F 5 , and Fq, which undergo both 
magnetic breakdown tunnelling and finite Bragg reflection (that is, 1 — P) are approximated 
to have the same g-factor with anisotropy This anisotropy is smaller than for orbits 
F 3 , which show only magnetic breakdown tunnelling without finite Bragg reflection, and which 
are approximated to have a common g-factor gp with anisotropy Since it is not possible 
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to identify unique gf-factor values given the multiplicity of frequencies, the values of g* and 
here represent parameters used for the simulation. 

The amplitudes of the F 2 and F 5 orbits have approximated to be equal, while experimentally, 
the amplitude of the F 2 orbit is slightly larger in contactless resistivity^^ and magnetic torque"^"^ 
measurements, probably due to additional effects such as small differences in the scattering 
time or effective mass. We note that the amplitude of the F 2 orbit is significantly larger in c-axis 
transport experiments.^^ 

We note that in Eqns. 2, 4, 5 the most general circular Fermi surface cross-section has been 
assumed, which is not specific to any Fermi surface reconstruction model. In order to extract 
the in-plane Fermi surface topology from the experimental data, a complete determination of 
the Fermi wavevectors would be required, for which future complementary experiments to 
the present quantum oscillation experiments are indicated. 

Additional angle-dependent damping factors considered. 

The inclusion of damping of the quantum oscillation amplitude from a finite lifetime of the 
quasiparticles would further not be expected to alter the conclusion of a staggered twofold Fermi 
surface geometry. This damping can be described in terms of a complex dispersion relation or 
effectively in terms of complex Fermi wavevectors, which are consistent with the symmetry of 
the full lattice potential both in their real and imaginary parts. Here the effect of quasiparticle 
scattering lifetime over the Fermi surface is included in terms of an imaginary component of 
the Fermi wavevector. The anisotropy of the amplitude arising from a complex Fermi surface 
(defined by Fermi wavevectors with real and imaginary parts) is expected to be the same as 
the original real Fermi surface. If the imaginary part of the warping were to dominate the 
real part then the amplitude damping factors would be given by expressions similar to Eqns. 3 
and 4 to leading order in warping, but with the outer Bessel function Jq replaced by the modified 
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Bessel function Jq, and AFq and AFi replaced by analogous parameters that measure the degree 
of damping. Our conclusion wherein the quantum oscillatory frequency spread is inconsistent 
with a fundamental neck and belly warping remains unchanged for the reasons that (1) for a real 
part of the warping sizeable enough to capture the frequency spread, a Yamaji resonance would 
occur, contrary to observation, whereas (2) if the real part were negligible then the amplitude 
would be governed by the very different function Jq, which is also contrary to experimental 
observation. 

Other amplitude damping factors that arise from random quenched sample inhomogeneities 
or magnetic field inhomogeneities would not be expected to show the symmetries of the under¬ 
lying lattice. Further damping factors such as additional damping within the vortex regime, or 
a damping factor of the form exp(—l/(cos 6')““^) where a is a variable parameter (after ref."^^) 
have been considered. However, such damping factors would yield an increasing amplitude 
suppression with increase in polar angle, and hence do not provide an explanation for the ob¬ 
served anti-resonance in quantum oscillation amplitude in the vicinity of a special polar angle 
^ 60 ° in underdoped YBa2Cu306.56- 

Negative Hall effect from nodal Fermi surface pocket from charge order. 

Of the various models proposed for the normal ground state of the underdoped copper ox¬ 
ides,examples of Fermi surface models that have been proposed in¬ 
clude An accurate Fermi surface model would need to explain 

the observation of quantum oscillations in underdoped YBa 2 Cu 306 +x on a background of a 
negative Hall coefficient (ref.^°), given that a positive Hall coefficient would be expected 
for the hole carrier doping in underdoped YBa 2 Cu 306 +x- In various previously proposed inter¬ 
pretations, the negative Hall effect is yielded by electron pockets at the antinodal locations of 
the Brillouin zone.^’ Our present work that locates the small Fermi pockets in the vicinity 
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of the nodal regions of the Brillouin zone, however, renders these interpretations unlikely. 

The nodal staggered twofold Fermi poeket in the new reeonstrueted body-eentred orthorhom- 
bie Brillouin zone that we propose to arise from staggered eharge order (sehematie in Extended 
Data Figs. 6 and 7; j-efs.^’^’^^’^°’^^), has opposite polarity to the original large Fermi surfaee in 
the starting primitive orthorhombie Brillouin zone in underdoped YBa2Cu306+x- In the high 
magnetie field limit, the Hall eonduetivity from a Fermi surfaee poeket — rather than being 
determined by the Fermi surfaee eurvature, as would be the ease in low magnetie fields — 
would instead be expeeted to be proportional to the number of states eontained within that 
poeket, with a sign that is negative for filled states and positive for empty states.A neg¬ 
ative Hall eoeffieient would therefore arise from the nodal staggered twofold Fermi surface 
from charge order,^’^’^^’^®’^^ providing an explanation for the previous puzzling observation of 
quantum oscillations in the negative Hall coefficient. 

Low quantum oscillation frequency. 

Very slow quantum oscillations are further seen over a large magnetic field range in contactless 
electrical resistivity measurements (Extended Data Fig. 3b). The separation between observed 
oscillation maxima (lower inset to Extended Data Fig. 3b) corresponds to a frequency of 90 ± 10 

T. 

A property of systems exhibiting magnetic breakdown is the occurrence of Stark quantum 
interference effects in electrical transport. Quasiparticles that encounter the magnetic 

breakdown junction (indicated in magenta in the upper inset of Extended Data Fig. 3) can take 
two possible paths (i.e. the inner orbit or the outer orbit), depending on whether or not magnetic 
breakdown tunnelling occurs. When the quasiparticles recombine at a second magnetic break¬ 
down junction they acquire a relative difference in Onsager phase (A^ = 27rAF/B), which is 
proportional to the difference in A;-space area between the two paths. Because this area is not 
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a closed orbit, it does not yield quantum oscillations in thermodynamic quantities such as the 
magnetisation or heat capacity — however, interference between the two paths can give rise to 
oscillations in electronic transport properties such as the electrical resistivity, or thermal con¬ 
ductivity. The simulation of the charge order model made using parameters in External Data 
Tables 1 and 2 yields AF ^ 90 T for the frequency corresponding to the small difference in 
area between the two magnetic breakdown junctions, suggesting the association of the observed 
low quantum oscillation frequency with the Stark quantum interference effect. Alternatively, 
this low frequency could correspond to the small hole pocket predicted in ref.^^ The existence 
of a small hole pocket in addition to the nodal electron pocket discussed in the main text could 
provide a possible explanation for the observed oscillatory component of the Hall resistivity.’’^ 

Anisotropy between in-plane and inter-plane effective hopping. 

The velocity anisotropy corresponding to a staggered twofold Fermi surface geometry may be 
estimated by means of a quasi-two dimensional energy band dispersion of a form consistent 
with the order of expansion in cylindrical harmonics in equation 4. In cylindrical coordinates 
(k, 0, kA and to an additive constant, this dispersion relation near the Fermi surface can be 
expressed in the form 

h^k"^ 

e =- h 2t* cos kzic 12) sin 20, 

2m* 

where the magnitude of the interlayer dispersion varies around the circumference. Here m* 
is an in-plane effective mass, kp is the effective Fermi wavevector, and t* is a c-axis effective 
hopping parameter. 

The ratio of the in-phase Fermi velocity to c-axis root mean square velocity v± is then 

v± ^AFt^oioid ^kpc 
where AFt^ofoid/Fo = 2t*/{h‘^kj/2m*). Taking Fq ^ 530 T, AFwoW ~ 15 T, 
kp = \j2eFo/h ^ 0.13 A~^ and c 11.7A, the velocity anisotropy ratio v\\/v± is of the order 


27 



of 100, a large value consistent with a nodal location of the Fermi surface pocket, given that the 
interlayer hopping is expected to be weakest in this region of the Brillouin zone.^^ Interestingly, 
the velocity anisotropy ratio is of a similar order of magnitude as the conductivity anisotropy 
ratio inferred for optical conductivity measurements in the pseudogap regime. 

The correlation length of the associated superlattice structureis required to be similar 
to or greater than the extent of motion of the cyclotron orbit in the interlayer direction, for our 
analysis arriving at the Fermi surface geometry in this work. The large anisotropy in velocity 
between in-plane and interlayer directions means that the cyclotron orbits in real space traverse 
an interlayer distance of only about 1 % of the in-plane cyclotron orbit radius, placing us in the 
required limit. 
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Extended Data 
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Extended Data Figure 1: Measured projection of the magnetic field along the crystalline 
c-axis of the sample. Circles indicate the maximum B cos 6 measured at 65 T (cyan) and 85 T 
(red), obtained by means of a projection coil, while the dashed lines represent fits to a cosine 
function. The angular error is < 0.2° for 9 < 66° and ~ 0.2° for 68° < 9 < 71°. 9 = 0°, 
1.3°, 11.3°, 12°, 16.3°, 18°, 21.3°, 26.3°, 31.3°, 36.3°, 38°, 41.3°, 45.2°, 46.3°, 48°, 49°, 49.4°, 
50.1°, 50.6°, 51.4°, 51.5°, 52°, 52.3°, 52.5°, 52.9°, 53.1°, 54.4°, 54.9°, 55.5°, 56°, 56.2°, - 
56.95°, 57.2°, -57.4°, -58.15°, 58.2°, -59.4°, 59.6°, 60.6°, 61.2°, -61.4°, 61.7°, 62.5°, 62.6°, 
-62.7°, -63.2°, 63.4°, 63.7°, -64.1°, 64.5°, 65.5°, 66°, 66.3°, 68.1°, 69.4° and 70.6°. Negative 9 
angles refer to a measured equivalent (180-101) angle as shown. 
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Extended Data Figure 2: Experimental quantum oscillations for different angles compared 
with simulations for a neck and belly model, a, Measured oseillations in the eontaetless 
resistivity, b, Simulated oseillations at the same angles and fields as a for two Fermi surfaee 
eylinders exhibiting a fundamental neek and belly warping, for parameters used in ref.^^ (listed 
in Extended Data Table 3) to simulate the restrieted experimental range within the dashed line. 
Data in a and simulations in b have been sealed by exp for visual elarity. c, Symbols 
represent the absolute value of the eross-eorrelation between the quantum oseillation data in 
a with a simple sinusoid exp {iF is matehed to the periodieity of the oseillations at 
9 = 38°, where a single frequeney dominates the measured quantum oseillations. Coloured 
lines indieate a simulation oc R„Rs for a staggered twofold model using parameters in Extended 
Data Table 1 (magenta), and a neek and belly model using parameters from ref.^® in Extended 
Data Table 3 (red). While the anti-resonanee ii?8ie vieinity of 9 = 60° yielded by the staggered 
twofold model is in good agreement with the experimental data, the striking Yamaji resonanee in 
the vieinity of 0 = 60° yielded by the neek and belly model is in marked eontrast to experiment. 
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Extended Data Figure 3: Quantum oscillations extending to 100 T magnetic field, a, 

Schematic illustration of the small and large Fermi surface pocket sizes (and quantum oscilla¬ 
tion Fourier frequencies) in underdoped YBa2Cu306+x and overdoped Tl2Ba2Cu06+5, respec¬ 
tively. b, Contactless electrical resistivity of YBa2Cu306.56 measured to 100 T, showing 
the resistive transition (at 20 T) and quantum oscillations. The dominant quantum oscilla¬ 
tions of 530 T can be seen to be superimposed on slowly varying oscillations (red line), which 
we extract in the lower inset by subtracting the dominant oscillations and a linear background. 
The slowly varying oscillations are consistent with a low frequency of 90 ± 10 T. The upper 
inset shows the bilayer-split pockets expected for charge order in which the difference in area 
between two magnetic breakdown junctions corresponds to a frequency of ~ 90 T. 
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Extended Data Figure 4: Contour representation of a staggered twofold Fermi surface 
model compared to the experimental data, a, Contour plot of the simulated quantum os¬ 
cillation amplitude for a staggered twofold Fermi surface geometry represented by equations 2 
and 4, using parameters in Extended Data Table 1 and shown in Fig. 4c. b, Contour plot of ex¬ 
perimentally measured quantum oscillation amplitude; good agreement is seen with the model 
in a. The quantum oscillation amplitude is indicated by the colour scale (in arbitrary units) in 
the reciprocal field-angle plane; for clarity the ordinate is given as tan Q. 


40 
















Extended Data Figure 5: Cross-correlation measured over different field ranges, compared 
to simulations for a staggered twofold Fermi surface geometry. Real eomponent of the 
eross-eorrelation between the quantum oseillation data over fixed ranges of B cos 9 for a range 
of measured 6 angles, with a simple sinusoid cos{2tiF/B cos 6* + 0). F and 0 are matched to 
the periodicity and phase of the oscillations at 0 = 38°, where a single frequency dominates 
the measured quantum oscillations. Black lines indicate the simulation (proportional to R^Rg), 
where AFtwofow ~ 15 T is the depth of the modulation, while square symbols indicate the 
experimental cross-correlation. 
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Extended Data Figure 6: Schematic of a nodal Fermi surface from charge order that stag¬ 
gered perpendicularly to the hilayers. a, Reconstruction of the Brillouin zone, with one 
instance of the pocket location indicated at the ‘T’ point in relation to the original Fermi surface 
(purple) and nodes in the superconducting wave function (from refs.^’^’^°’^'). Here the concen¬ 
tric arrangement of Fermi surfaces arises from bilayer splitting. The in-plane shape of the Fermi 
pocket shown here is an illustration based on the model in refs.;^’^’^°’^' the corrugation shown is 
significantly accentuated for clarity, h, A three-dimensional view of a to illustrate the staggered 
twofold geometry. While this schematic assumes achirality, a chiral model is not ruled out, as 
for instance proposed in ref.,^^ where the form of order breaks mirror symmetry within each 
plane. 
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Extended Data Figure 7: Schematic of the Brillouin zone cross-section, showing magnetic 
breakdown orhits in a charge ordering scheme. A cut through the = 0 plane of the 
Brillouin zone shows the six possible orbits resulting from magnetic breakdown tunnelling in a 
bilayer charge ordering scheme/^’an illustrative Fermi pocket shape similar to Extended Data 
Fig. 6 is shown, a The two Fermi surface cross-sections of frequency Fi = — 2AFspiit and 

Fq = Fq + 2 AFspiit that can result from bilayer splitting with in-plane ordering wave vectors 
qI = 27r(±^, 0) and Q| = 27r(0, ±y). The T and T symmetry points of the body-centered 
orthorhombic Brillouin zone of the charge order superstructure are depicted in blue. The gap 
separating bonding and antibonding surfaces is expected to be smallest at the nodes.Panels h, 
c, d and e show the range of possible magnetic breakdown orbits, F 2 = Fq — AFgput, F 3 = Fq, 
F 4 = Fq and F^ = Fq + AFsput, as listed in Table 2. 
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parameter 

description 

value 

Fo 

quantum oscillation frequency 

534 T 

^.^twofold 

staggered twofold warping frequency 

15 T 

AFsplit 

bilayer splitting frequency 

90 T 

mJi 

quasiparticle effective mass 

1.6 TTie (fixed) 

Bo 

magnetic breakdown field 

2.7 T 

A 

Dingle damping parameter 

132 T 

Oq 

absolute amplitude 

0.131 MHz 

9\\n 

g-factor 1 

2.12 


g-factor 2 

0.35 


g-factor anisotropy 1 

1.43 

^0 

g-factor anisotropy 2 

0.15 


Extended Data Table 1: Model parameters for a staggered twofold Fermi surface model. 

Tabulated values used in equations 2 and 4 to simulate the quantum oscillation waveform, yield¬ 
ing excellent agreement with experiment as a function of B, 9 and 0 (Fig. 4 and Extended Data 
Figs. 2c and 4). The effective mass is taken to be a fixed quantity, having been determined in¬ 
dependently from temperature-dependent measurements."*^ The parameters are the same for all 
the orbits, except for those denoted by subscripts □ and <0, which each correspond to a subset of 
orbits as defined in the text. Because of the multiple frequencies in the model, it is not possible 
to uniquely identify the g-factor; the values of g* and here represent parameters used for the 
simulation. 
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orbit Fj 

A, 

-Rmb 

Fi = Fq — 2APspiit 

1 

1 

to 

F 2 = Fq — APgpiit 

4 

-P(l-P) 

F 3 = Fq 

2 

p2 

P 4 = Fq 

4 

-P(l - P) 

F 5 = Fq + APgpiit 

4 

-P{1 - P) 

Fq = Fq + 2APspiit 

1 

{i-py 


Extended Data Table 2: Magnetic breakdown amplitude damping for the orbits shown in 
Extended Data Fig. 7. The magnetic breakdown network corresponds to a split Fermi surface 
geometry from charge order, as shown in Extended Data Fig. 7. A high magnetic breakdown 
tunneling probability (P) causes the amplitude of the Pi and Fq orbits to become very weak, as 
seen in Fig. 3a (ref.^^). 


parameter 

description 

value 

Fi 

quantum oscillation frequency 1 

478 T 

F 2 

quantum oscillation frequency 2 

526 T 

^.fjieck—belly, 1 

Fermi surface corrugation frequency 1 

37.7 T 

'^-^neck—belly,2 

Fermi surface corrugation frequency 2 

3.5 T 

7i 

quantum oscillation phase 1 

3.5 

72 

quantum oscillation phase 2 

1.1 

ml 

quasiparticle effective mass 1 

1.5 TTle 

ml 

quasiparticle effective mass 1 

1.7 me 

Ai 

Dingle damping parameter 1 

127.8 T 

A 2 

Dingle damping parameter 2 

159.8 T 

ai 

absolute amplitude 1 

0.238 MHz 

0-2 

absolute amplitude 2 

0.339 MHz 

9h 

g-factor 1 

1.40 

9\\2 

g-factor 2 

1.88 


Extended Data Table 3: Model parameters for a fundamental neck and belly Fermi surface 

model. Tabulated values from ref.^^ used in equations 1 and 3 to simulate the waveform of 
the oscillations in Fig. 2c and Extended Data Fig. 2b and c for two cylinders with isotropic g- 
factors. The absolute amplitudes ai and 02 listed here have been adjusted to match the amplitude 
in the present experiment. 
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